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Introduction 


This presentation will go into detail concerning one specific application: 
moving body simulations of the MSL capsule for the purpose of ballistic 
range test comparisons and MEDLI 2 design. 


¢ MEDLI: Mars Science Laboratory . ——- a ; 
Entry, Descent, and Landing Instrument oo _ JE® 


¢ MEDLI 1 was for the MSL mission. . 
¢ MEDLI 2 is for the Mars 2020 mission. i 


~ 


fa) OML geometry 


Simulations were performed using the OVERFLOW Navier-Stokes solver 
in support of the 90mm ballistic range test. The test was conducted 
earlier in 2016, but the data are still being reduced. This presentation 
will cover the pre-test CFD that was performed. 


For dynamic simulations with a pitching capsule, overset grids were an 
obvious choice to enable motion of the capsule without affecting the 
resolution in the shock and wake. 
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Grid Generation 


To support initial static effort and future dynamic and 6-DOF runs, a grid 
script system was made to enabled grid motion and shock adaptation. 


¢ Chimera Grid Tools (CGT) allow for repeatable, dynamic generation of surface 
and volume grids. 


¢ Grid spacings initially based on Orion MPCV capsule grid scripts. Incorporating a 
decade of grid refinement studies for capsule shapes. 


¢ Surface grid discretization was further refined to more closely match the legacy 
LAURA grids used by the MEDLI team. 


Shock Box 


Wake Box 


Automatic 
Off-Body Grids 
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Grid Generation 


OVERFLOW’s build-in Domain Connectivity Function (DCF) used to 
perform hole cutting and interpolation for static and dynamic cases. 


¢ DCF enables arbitrary orientation of capsule prior to simulation and enables moving 
body. It also enables automatic off-body grid generation. 


¢ Grid is rotated about the capsule CG to vary attitude. 


a = 0° a = 30° 


Rotation 
——» 
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Grid Generation 


OVERFLOW’s build-in Domain Connectivity Function (DCF) used to 
perform hole cutting and interpolation for static and dynamic cases. 


¢ DCF enables arbitrary orientation of capsule prior to simulation and enables moving 
body. It also enables automatic off-body grid generation. 


¢ By using DCF, we were also able to use OVERFLOW’'s grid adaption to investigate 
grid sensitivity and shock resolution. 


a = 0° a = 30° 


Rotation 
Adaptation 
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Effects of Adaptation 


Simple refinement study showed diminishing returns with grid adaptation. 


¢ Only affect is on heatshield pressure. 
¢ Base pressures are constant with additional refinement. 


CPU 
Mach 3.0 Relative 
a=30° ma Cn Cm,x=0 helt Cost 
submission) 
2 Levels 1.2739 - -0.0917 - 0.0756 - 177 168% 


Adaptation 


Love 1.2723 0.13% | -0.0918 0.01% | 0.0756 0.03% 147 139% 
Adaptation 


No- 1.2697 0.33% | -0.0914 0.40% | 0.0753 0.37% 106 100% 
Adaptation 


One level of adaption was used for the static simulations. No 

adaptation was used for the dynamic simulations. 

¢ Not used for dynamic simulations in order to minimize cost since benefit was 
small. 
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OVERFLOW Inputs 


Simulations use OVERFLOW version 2.2k. 

¢ Turbulent, viscous Navier-Stokes solver 

¢ Using SST-DES turbulence model (NQT=205, IDES=1) 

¢ Using HLLE++ numerical fluxes (IRHS=6, DELTA=5.0) 

¢ Implicit SSOR time advancement (ILHS=6) 

¢ Adaptation on undivided differences (ETYPE=0, MAX_SIZE=80.e6) 
¢ Time-accurate, unsteady flow (DTPHYS=0.04, NITNWT=4-8) 


- Sensitivity study was performed to determine the timestep. 
- Cases show 3-4 orders of magnitude drop in average residual during sub-iterations. 


¢ Viscous, Isothermal wall temperature of 540 °R 


At JSC, we have developed a capable OVERFLOW run manager, overlst’. 


¢ It is a collection of Perl modules that interface with CGT and OVERFLOW as well as 
several compute clusters. 


¢ Relevant portions have been augmented to handle DCF and moving body simulations. 
¢ This run manager handled the nearly 100 preliminary, static, and dynamic simulations 
used for this work. 
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Cas [-] 


Results - Static Aerodynamics 


Static aerodynamics were integrated from static simulations: 
¢ Sea level conditions (simulating ballistic range) 
¢ Mach numbers considered: 1.5, 2.0, 2.5, 3.0 °Q°<a< 30° 


Goal was to develop static aerodynamic curves to compare to existing 
models and the future dynamic work. 


¢ Mach 1.5 reduction in axial force at low a corroborated by other sources. 


Mach 3.0 -®@ Mach2.5 -®© Mach2.0 -® Mach1.5 -& 


1,05 0.1 
0.09 
1.5 0.08 
1.45 oe 
—~ 0.06 a 
1.4 S. 005 2 
O 0.04 O 
1.35 0.03 
13 0.02 
0.01 
1.25 | | | | | 0 
O 5 10 15 20 25 30 O 5 10 15 20 25 30 O 5 10 15 20 25 30 
o [°] o [°] o [°] 
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1-DoF simulations were performed looking at pitch damping from an 
initially static MSL at varying angle of attack. 
¢ Initial amplitudes of 5°, 10°, 15°, and 30° at Mach 1.5 and 3.0 


¢ Same conditions and numerics used previously, but no adaptation. 
Mach contours from case initially at 30°. 
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Aerodynamic Model 


Assume linear response in dynamic derivatives with respect to pitch- 


rate (q) and that dynamic derivatives are only a function of a and not q. 


¢ Based on the data, there is evidence of non-linearity in dynamic derivatives. 


Co = Cg(a) + Co,q(@) xq 
SKS s/ 
static dynamic 


Based on the above, the dynamic simulations can be used to get 
estimates of static and dynamic aerodynamic coefficients. 


¢ Least-squares fit through available data at angles of attack. 
¢ Details are shown on following slides. 


Run strategy begins with converged static simulation and then allows 
motion. 


¢ First half period of oscillation thrown out to eliminate effect of transients. 
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Method of Calculation 


a— a 0° =10.0° a=25.0° 8 a=10.0 © 
Cn — 30° =5, = az5 a=20, 
0.04 
20° 
0.02 
10° 
0 oe & r 
5 oO 
oO 
-0.02 “10 
-20° 
-0.04 
-30° -30 
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 
Time [s] Time [s] 


From time history of aerodynamics and angle of attack, identify all 
points at a discrete angle of attack. 


¢ Increments of 0.5° were used to discretize attitude. 
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-30 
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 
Time [s] 


a [°] 


Method of Calculation 


C., [-] 


-0.02 


P @ «a=25.0° @ a=-100° © a=-25.0° © 
é 2-5, as { 


0.04 


-0.004 -0.002 0 0.002 0.004 
Pitch Rate (qD)/(2V) [-] 


For visualization, these data can be replotted versus pitch rate. 


¢ Inherent symmetry in the geometry and problem allow for mirroring all data at 
negative a to increase the effective number of points at positive a. 
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Method of Calculation 


-0.004 -0.002 0 0.002 0.004 -0.004 -0.002 0 0.002 0.004 
Pitch Rate (qD)/(2V) [-] Pitch Rate (qD)/(2V) [-] 


For visualization, these data can be replotted versus pitch rate. 
¢ Remember the model: C¢ = Cy (a) = C'y.q(a@) x qd 
——_—F i 


b max 


Slope of resulting least-squares fit is Cg 7(a) ; intercept is C's (@). 
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Results - Dynamic Simulations 


Varying the initial pitch amplitude affects the pitch rates that are 


observed at a given a. If Cm is linear as the model assumes and only a 
function of a and extracted aerodynamics from all four simulations 


should agree. 1-DoF ag = 30° 1-DoF ag = 10° 
1-DoFag=15° — 1-DoFag= 5° — 


30 


20 


a [°] 
© 


-30 
-0.004 -0.002 0 0.002 0.004 


Pitch Rate (qD)/(2V) [-] 
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Results - Dynamic Simulations 


To compare the validity of the method, we compare integrated 
aerodynamics from the non-moving simulations to those calculated 
with the curve fit of the moving simulations. 


Mach 3.0: 
¢ Fit with the 30° data compares very well to the non-moving results from earlier. 


¢ Fit with other angles is comparable, but noticeable deviations in Cas (non-linear?). 
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Results - Dynamic Simulations 


NES eee 
Sent e Landi vg 


Mach 1.5: 


¢ Only three cycles of motion versus six at Mach 3.0. Lower speed, generally more 
difficult and less stable. 5° simulation is noticeably undamped. 


¢ Data from moving simulations not as smooth as for Mach 3.0. Ca; shows larger 
deviations. For Cn,s and Cm,s the agreement is still good between static and 
dynamic simulations. 


Static @ 1-DoF dg = 30° — 1-DoF a = 10° — 
1-DoF dg = 15° — 1-DoF a = 5° — 
0.005 


Cas [-] 
CNn.s [-] 
Cm,s [1] 

o 

Oo 
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Results - Dynamic Simulations 


A ballistic range test was preformed for the MSL capsule ~10 years ago. 


¢ It was performed at a different scale as the current work (65mm versus 90mm) 
and using slightly different mass properties. 


¢ Ballistic range reductions for dynamic terms are quite different than the linear 
model and curve-fit that was shown previously. Portions of several runs are 
combined and data are fit using multiple coefficients. 


Cm, = 6) Sg Clase 9 Crag ale Cong m (Mi _ Mref) =F C 


Mgq,m2 


V2+V? 
€ = sin(ar) = 72 Meret S25 


¢ Regardless, our hope was that this data could give us an order of magnitude 
comparison to the CFD. 


¢ Data taken from AIAA 2009-3917: “Dynamic Stability Testing of the Mars Science 
Laboratory Entry Capsule” 


(M; _ Mrep)? 
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A ballistic range test includes placing 
a scale model inside of a sabot ata 
desired initial orientation. The sabot 
and model are place into a gun and 
fired down an instrumented range. 


¢ Onboard instrumentation 
¢ Photogrammetry stations 
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Results - Dynamic Simulations 


Order of magnitude is similar. Some observations: 


¢ Variation in the results from the simulations imply that our model is not perfect or 
that Cm,q is also a function of q. 


¢ At Mach 3.0, the data shows good agreement with the range data using the data 
from the 10° and 15° simulations. 


¢ Mach 1.5 show a reduction in damping as compared to Mach 3.0. Do not agree 
as well with the range data. 


Ballistic Range - - 1-DoF a = 30° — 1-DoF a = 10° — 

1-DoF a = 15° — 1-DoF dg = 5° — 
0.2 { 
0.1 0.8 
0 0.6 
Zo ——-_ Mach 3.0 — 0.4 
2 = 0.2 
gs -0.2 ; 
03 Mach 1.5 ————> Qs 
-0.4 -0.4 
-0.5 -0.6 

0 5 10 15 20 25 30 O 5 10 15 20 25 30 
a [7] a [°] 


13th Symposium on Overset Composite Grids and Solution Technology 21 


Conclusion and Future Work 


Ability to couple CGT grid scripts and OVERFLOW DCF was crucial to 
enable rapid development of grids and solutions. 


Method with rotating body-fitted, viscous grids showed agreement with 
static data and similarity with previous ballistic range data. 


Entire process (static, dynamic, reduction, and comparison) performed 
with non-overset solver (US3D) and agreement was very good. 


¢ More details can be found in conference paper: 
‘Dynamic CFD Simulations of the MEADS II Ballistic Range Test Model’ (AIAA 2016-3243) 


2016 ballistic range data is currently being reduced for additional runs. 
¢« New cases being run to match the range conditions. 


¢ Plan is to perform 1-DoF prescribed and 6-DoF free simulations for comparison using 
OVERFLOW. 


¢ Data reduction strategy will be augmented to include method used in the ballistic 
range tests. 
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Backup 
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Moving Body, 6-DOF Inputs 


prescribed motion input. 


Mass Properties for model: 

¢ Ix = 916460.86 [g*mm?] 

¢ ly = 619454.59 [g*mm?] 

¢ Iz = 656597.73 [g*mm*?] 

¢ Mass = 1369.76 [g] 

¢ (X,Y,Z)cc = (-0.3002D, 0.0, 0.0) 
Simplified for initial work to: 
¢ Ix = 916000 [g*mm7] 

¢ ly = 638000 [g*mm7] 

¢ Iz = 638000 [g*mm7?] 

¢ Mass = 1370 [g] 
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Grids system and inputs for OVERFLOW updated to allow for 6-DOF or 


